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ABSTRACT 

I present a general family of dynamical models with simple analytical 
potential-density pairs suited to model galactic bulges and nuclei with double 
power-law radial density profiles and an optional central black hole. Analytical 
expressions for the potential and velocity dispersion of these models are given 
in the spherical case. Many previously known analytical spherical models, 
including also the recent 7/?y-models by Dehnen (1993) and Tremaine et al. 
(1994), are special cases of this family. This family also forms a complete 
set for constructing general galaxy models or solving Poisson's equation in 
the non-spherical case. In particular, a generalized Clutton-Brock (1973) and 
Hernquist-Ostriker (1992) orthogonal basis set is given. 

Subject headings: celestial mechanics, stellar dynamics -galaxies: kinematics and 
dynamics - galaxies: nuclei, methods: analytical 

1. Introduction 

Recent high resolution photometric data from the Hubble Space Telescope (HST) 
(Crane and Stiavelli 1993, Kormendy et al. 1994, Ferrarese et al. 1994) reveal three basic 
features in the mass distribution of galactic nuclei: (1) a central cusp in radial density profile, 
p oc with < 7 < 3; (2) a markedly non-spherical shape (fiattened/triaxial/lopsided) in 
density distribution (except for perhaps the nearly spherical M87 nucleus); (3) a possible 
massive central black hole, which is suggested by the cusped light profile (Young 1976). 
Of these systems, M32 and M31 are two examples. The former has a fiattened E3 nucleus 
with a central cusp 7 ~ 1.5. The latter has a double-nucleus, where the fainter peak (P2) 
appears to be cusped in the high resolution HST observation and the off-centered brighter 
peak (PI) may generate a significant lopsided perturbation (m=l mode) in the potential 
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(Lauer et al. 1993). For both a massive central black hole of roughly a few times 10^ Af© 
and 10'' Af© respectively has been proposed based on in detailed axisymmetric models of 
kinematic data (Dehnen 1994, Qian et al. 1995); the results are in qualitative agreement 
with earlier spherical models (Dressier and Richstone 1988). 

Modelling the structure of nearby galactic nuclei can help us constrain their formation 
history and understand their relation with distant active galactic nuclei. A new class of 
models are needed to account for the above three basic features in recent observations. 
The simple parametrization first proposed by Hernquist (1990) appears to be an attractive 
choice. It is a dimensionless spherical volume density model with three free parameters 



where C is a normalization constant. This model parametrizes the volume density as a 
general double power-law with slope —7 at radii much smaller or —fi at radii much bigger 
than the break radius r = 1. The third parameter, a, is a measure of the width of the 
transition region: the bigger a is, the wider the transition region. Similar parametrizations 
have been used to fit observations of galactic nuclei (Kormendy et al. 1994, Crane and 
StiaveUi 1993). 

Hernquist (1990) proposed the (a, /3, 7)-models only as a density model. Their 
dynamical properties are not known analytically except for the Hernquist bulge and a few 
other well-known special cases (see Table I), and, more recently, a narrow sequence due to 
Dehnen (1993) and Tremaine et al. (1994). The latter authors found that the spherical 
(1,4,7) subset has many good analytical properties, including analytical potential, intrinsic 
velocity dispersion for all real values of 7 between and 3, and for 7 = 0, 1,2 analytical 
projected light and dispersion as well, all of which remain analytical when a central black 
hole is introduced. I shall later on call this (1,4,7) subset as 7/?y-models, which relate to 
the ?y-models by Tremaine et al. (1994) by 7 = 3 — 77. Note throughout the paper, 7 always 
denotes the inner power-law slope in the volume density, fi the outer power-law slope. All 
models have a finite mass. The total luminous mass and the gravitational constant G are 
always set to unity. 

In this paper, I show that the potential of the whole set of (a, /3, 7) density models can 
be computed analytically, and the isotropic velocity dispersion of most of these spherical 
models with or without a black hole can be given analytically. Table I shows that the family 
includes many well-known spherical models, e.g., the Hernquist bulge (Hernquist 1990), the 
Jaffe model (Jaffe 1983), the Plummer model (1911), the modified Hubble profile (Binney 
and Tremaine 1987), the modified isothermal dark halo model (Sackett and Sparke 1990), 
the perfect sphere model (de Zeeuw 1985a, b) and the recent 7/77 series (Dehnen 1993, 



(«,/3,7), 
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Tremaine et al. 1994). But most of the models in the three parameter spherical family 
are previously unknown. Some have very simple potentials, like the a model and the /3 
model. The (n,3 + ^,7) three paramemter family also have analytical velocity dispersions. 
Besides a series of models included in the known 7/77 set, there are also two new models 
with analytical distribution functions. 

The simple analytical properties of the spherical models make them useful for i) 
deriving the mass-to-light ratio (Tremaine et al. 1994) as traditionally done by core fitting 
with a King model, ii) setting up equilibrium conditions for N-body simulations to examine 
the secular evolution of galactic nuclei (Quinlan et al. 1995) and iii) testing numerical 
models (Dehnen and Gerhard 1994). While modelling the light and kinematic data of 
galactic nuclei is generally a numerical problem, analytical potentials can still serve as a 
solver for Poisson's equation and speed up the force calculation in N-body simulations 
(Hernquist and Ostriker 1992, Hernquist et al. 1995) and in building orbit models (Zhao 
1994). 

Section 2 shows that the potential of the general (a,/3,7) density models can generally 
be written analytically for both the spherical and the non-spherical cases. Section 3 presents 
the exact analytical expressions and asymptotic expressions for the dynamical quantities in 
the new spherical models. A brief summary is in Section 4. Many less important details of 
the models are in the Appendix. 

2. Analytical potentials of (a, /3, 7)-models 

Table I lists various new analytical models and known models which are subsets of the 
spherical (a, /3, 7)-models. One can see that the 7/?y-models belong to a three paramamter 
(n,3 -|- ^,7) family which has analytical dispersions; n and k are any natural numbers. The 
latter is one of the four spherical families whose potentials can be expressed in terms of 
elementary functions. The potential can always be reduced to the incomplete Beta-function 
for both the spherical and the non-spherical case if the radial profile of each angular term 
in the density model is a double power-law with generally different slopes. 
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2.1. Two simple examples 



Before showing the details of the general models, I give two simple examples. It is 
straightforward to apply the Poisson's equation to the following two potentials, 

^(r) = --(l ^ 



r' (1 + r)/^-3" ' ' 

and ^ 

and to show that their corresponding densities have double power-law profiles. The former 
makes up a new analytical potential-density pair which I call the (1, /3, l)-models, or in 
short the /3-models. These models have a fixed inner density profile of but a variable 
outer slope /3 in contrast to the 7/77 models, which have a fixed outer profile r""* but a 
variable inner profile. The latter analytical potential-density pair is the (a, 3 -|- ^,2 — ^) 
models, which I call the a-models. They fill the gap between the familiar Plummer model 
(Plummer 1911; Binney & Tremaine 1987) with a = | and the Hernquist model with a = 1. 

The limitations of the a-models, /3-models and the 7/?y-models are obvious when it 
comes to fitting observations. Generally one cannot expect a good fit for all three regions of 
a spherical nucleus, namely, a central cusp, a power-law tail and the transition region, with 
only one free parameter. The central bulge and nucleus region of our Galaxy, for example, 
has a double power-law of slope —1.8 and —3.7 inside and outside 1 kpc respectively 
(Sellwood and Sanders 1988), which cannot be well-approximated by any of these simple 
models. And like the central regions of M32 and M31, the system is far from being spherical 
based on recent near infrared maps of the Galactic plane from COBE (de Zeeuw 1993, 
Dwek et al. 1995). For these reasons, I will focus on the more interesting subsets of the 
spherical (a, /3, 7)-models which have simple analytical expressions for dynamical quantities 
but fiexible density profiles. The generalization to the non-spherical case is straightforward. 



2.2. General analytical spherical subsets 



Consider the potential for the spherical (a,/3,7) density models. Substituting 
equation 1 to equation (2-122) of Binney and Tremaine (1987), one finds that 

r-"'-^(l + r^/")-('^-"')"Jr + r-^y r2-"'(l + r^/«)-('^-"')«Jr) . (4) 
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After a change of variable from r to x, 

Y = ^ = ^^ (5) 

the integrals in the expression for the potential can be written in terms of the incomplete 
Beta-function, 

$(r) = -4:7rp{r)fo,o{r)r\ (6) 
r /X ^ aB{co - go,go,x) aB{co -po,po,l - x) 

where the constants po, <?o and Cq are certain combinations of a, /3 and 7 given in the 
Appendix. From the properties of the incomplete Beta- function (see the Appendix), one 
can prove that the model potential $ reduces further to elementary functions if ONE of the 
following four pairs of numbers is a pair of natural numbers: 

(co - po, qo), [po, Co - qo), [po, qo) and (cq - po, Cq - qo). 

Among the four spherical families, the 7/?y-models, /3-models and a-models shown 
earlier belong to the first three families with the natural number pair being (2, 1), (1,2) 
and (1,1) respectively. Other explicit expressions will be given only for the family with 
(co — Po,<?o) being a pair of natural numbers (n -|- A;, A;), which is the most interesting 
family because it has analytical velocity dispersion as well as potential. See equation 17 of 
Section 3.1.. 



2.3. Generalization in shape and radial profile 



The analytical expression for the potential is also readily generalizable to the non- 
spherical case. This requires that the density of the observed system in polar coordinates 
(r, ^, (f>) can be decomposed into spherical harmonics y;^m(^,V') 

pir, ^, V) = E V^n™(^, i^)pUr), (8) 
and each angular component (l,m) has a double power-law radial profile 

where [ai^fii^^i) are three parameters to be adjusted to fit observation at the transition, 
inner and outer regions respectively. The corresponding potential has an analytical 
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expression very similar to that of the spherical model, 



$(r, 9, tjj) 




l,m 



(10) 



fi,m{r) 



(11) 



In the spherical case, equation 8, 10 and 11 become equation 1, 6 and 7, as the coefficients 
Ci^m are all zero except for a normalization constant Co,o = C*, which is given in the 
Appendix. The subscript / or (/, m) is sometimes omitted in the spherical models, where 
only the (/,m) = (0,0) term is relavent. To simplify notation, a subscript m is also omitted 
for the parameters p;, qi and q, which can depend on the azimuthal number m. The 
potential again reduces to elementary functions if one of the following is a pair of natural 
numbers 



For an observed system whose the density profile can not be well-fitted by a single 
double power-law (equation 1), it is still of numerical interests to fit the profile with a 
few double power-law components because the total potential will then be written as a 
linear combination of several explicitly known analytical components. The computation of 
force can be faster than fully numerical methods, particularly if the observed system has a 
double power-law profile to the zeroth order. The same is true for the angular expansion 
of the density distribution of a general non-spherical nucleus. In fact the whole family of 
(a, /3, 7)-models and their spherical harmonics terms form a complete basis set for solving 
the Poisson's equation. The orthogonal basis functions shown next are linear combinations 
of the subset of the (a, /3, 7)-models with (pi^qi) being natural numbers. Although the 
whole (a,/3,7) set is degenerate and non-orthogonal, it can still be useful in new methods 
to solve Poisson's equation with a general non-orthogonal basis set (Saha 1993). However, I 
will focus on the orthogonal set, which is easier to apply. 

2.4. a-model and orthogonal basis functions for Poisson's equation 

An efficient Poisson solver is crucial in massive N-body simulations since the 
accelerations of particles need to be computed at each time step. One of the most efficient 
methods is the Self- Consistent Field (SCF) method by Hernquist and Ostriker (1992). In 
this method, the potential and the density are written as a superposition of orthogonal 
basis functions, which are solutions of the Poisson's equation. The gain in speed depends 
on carefully choosing the lowest order term, so that one is able to resolve the density 



(q - pi, qi), {pi, ci - qi), {pi, qi) and (q - pi, q - qi). 
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and the potential in the least number of expansion terms. The lowest order term in the 
Clutton-Brock (1973) expansion is the Plummer model, which is suitable for systems with 
a finite core. Hernquist and Ostriker (1992) generalize the method so that the lowest order 
term can be the Hernquist model, which well approximates the i?4 law distribution of 
ellipticals and bulges. Since galactic nuclei neither have cores nor follow a unique Ri law, 
these two known expansions are clearly not always the optimal. 

One can easily build an optimal basis set by choosing the lowest order term to be the 
more general a-model shown in the Section 2.1.. Depending on whether the interesting 
region is the central cusp or the halo, one can choose the value of a according to the inner 
power-law slope 7 = (2 — ^) or the outer slope /3 = (3 + ^). The orthogonal basis functions 
are given as follows in notations similar to Hernquist and Ostriker (1992). 

( a i\ Kn,iV^Yi,m{e,i^) r' 

<5n,iAr, 0, V-) = (^i^^i/.yi. ^n (0 , (13) 

where G'^{^) is the Gegenbauer polynomials, ^ is given in equation 5 and the constants 
Kn,h ^n,i and w are given in the Appendix. One recovers the Hernquist-Ostriker expansion 
and the Clutton-Brock expansion by letting a = 1 and | respectively. Such expansions 
have been applied to approximate the potential of the boxy Galactic bar (Zhao 1994) based 
on the COBE observation (Dwek et al. 1995). The application of such models in the fully 
parallel Self- Consistent Field code (Hernquist et al. 1995) remains to be seen. 



3. Dynamics of analytical spherical (a, /3, 7)-models 



A new family of dynamical models that are analytical in both velocity dispersion 
and potential with or without a central point mass can be obtained within the spherical 
(a, /3, 7)-models. This property allows us to study the dynamics of a cusped nucleus under 
the infiuence of a central black hole analytically under the assumption that the phase space 
density depends only on the energy integral so that the velocity dispersion (T^(r) of the 
system is everywhere isotropic. Much of the calculation can be done in the same fashion as 
in Tremaine et al. (1994). 
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3.1. The analytical (n,3 + ^,7) subset with/without a black hole 

I find that in many cases cr(r), $(r), M[r) and p[r) can all be written as elementary 
functions of r. The most general and useful subset of these analytical models is the 
(n,3 + ^,7)-models with an optional central point mass, where n and k are any natural 
numbers and 7 is a real number between and 3. Also any linear combination of these 
models (with the same value of a) is also analytical. These models are less restrictive than 
the (1,4,7) models, and greatly broaden the range of analytical self-consistent models with 
black holes. The projected density and dispersion are also analytical if one restricts oneself 
to the (1, 3 + A;, 0), (1, 3 + A;, 1) and (1, 3 + A;, 2)-models. See Table 1. 

For the (n, 3 + ^, 7)-models with a constant mass-light ratio T and a central black hole 
of mass niBH} 

p{r) = Cx-""'(1 - X)"^ (14) 
M{r) = mB/j + 47raC5(a(3-7),a(/3-3),x), (15) 

a{fj-3)-l 

= mBH+ E «»-x"^'""'^+\ (16) 

8 = 

a(f3-2)-2 

Hr) = J: fe8^„(2-.)+8(x), (17) 



\r) = ^Y^^/^ ^XiM(r,(xi))xr^^^"^-^(l-Xi)"('^+^)-\ (18) 



(7 r 



^ 2a(f3-l)-2 a(f3+l)-l 



p{r) 



d^S2a{i-f)+z{x) + rriBH E ej5'„(_i_^)+j(x)), (19) 



8 = 8 = 



1 



"(2-7) 



m = ml-2H.iR)^ (20) 

8 = 



= t7^('E' 9^V^t2.)JR) + ^BHT.h.vtUR)), (21) 

V / 8 = 8 = 

where the coefficients a^, bi, di — hi and the functions Si{x), U-"'^ and V^-^"^ are given in 
the Appendix, x is given in equation 5. The above formulae can be derived by changing 
variable from r to x followed by integrating a polynomial expansion with respect to x- For 
the (1, 4, 7)-model, the above reduces to the formulae in Tremaine et al. (1994). 

There are many simple and plausible families in the analytical (n, 3 -|- ^, 7)-models, for 
example, the (1, 4, 7)-model, (1, 5, 7)-model, (1, 6, 7)-model, (2, |, 7)-model, (2, |, 7)-model, 
(2, Y7 7)-niodel. As an illustration, I give the closed form of these analytical properties for 



a(2f3-3)-l af3 
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the (1,5,7) models, which have a steeper outer power-law than the 7/77 model. 

Pir) - ^^-^^^^-'^ (22) 
^ - 47rr"'(l + r)5-7' ^^^^ 

M{r) = mB// + j/'-"'((4-7)-(3-7)j/), (23) 

m = -(4 -7)^2-7(2/) + (3 -7)^3-.(j/), (24) 

r 

^'(0 = --^-^[(4-7)5'2-27(j/)-(23 - 67)5'3-27(j/) + 5(10-37)5'4-27(j/) 

(1 - y) 

-10(7 - 2^)Ss-2-y{y) + 5(10 - 37)5'6-27(j/) - (19 - Q^)Sr-2-y{y) + (3 - ^)Ss-2-y{y) 

+mBH{S-i-y{y) - 5S.^{y) + 10S^.^{y) - 10S2-y{y) + 5S3-y{y) - S,.^{y))i (25) 

where J/ = X = r/[l-\-r). Although the expression for cr^ is long, it involves no cancellation of 
large terms and hence is stable when evaluated numerically; each Si{y) term is well-behaved 
mathematically and is finite except at the origin y = r = for z < 0. 

Analytical phase space densities are also available for somewhat restrictive but plausible 
models without a central black hole. These are related to the analytical potential-density 
pairs through the Eddington's formula (Binney and Tremaine 1987). Besides the Plummer 
model and the (1, 4, 2 — ^) series for all natural number values of n (Dehnen 1993, Tremaine 
et al. 1994), I find two more models with this property: the (1,5, l)-model with 

9 ^ yi(33 + 70e-114e2 + 32e3) 

J \^) 



2 



V2Tr^' 3(l+4e)(2-e; 

3 , 3'\/c \i / \ 

-|- arctan 2e H — — (arctan , -|- arctan , )J, (26) 

2(2 - e)^/^ v2 - e v2 - e 



and the (2, -, -)-model with 



27 y^(18 - 69e - 30e2 - 40e3 + 16e4) 



64y27r2' 97r(l-e)4 

6 — 27e -|- 56e^ , 2arcsinA/e,, 
+ 6(1 -eF' " + <2'' 

where — e is the energy integral. As a nucleus with an adiabaticly grown central black hole 
would have a central cusp of power | (Binney and Tremaine 1987), the (2, |, |) model can 
be of some interest. The (1,5, l)-model is very similar to the Hernquist model and can be 
useful for bulges with density fall-off steeper than the R^^'^ law at large radii. For alternative 
approaches to building models with analytical distribution functions, see Dejonghe (1984). 

It can be easily shown that all of the (a, /3, 7)-models have a positive definite 
distribution function f{E) > based on the Eddington's formula (Binney & Tremaine 1987) 
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since > 0. The stability ol these models is unclear. A sufficient condition for stability is 
'^'^j^^ < or more directly < (e.g., Binney and Tremaine 1987). The models are likely 
to be stable in the absence of a black hole: one can show that the sufficient condition is met 
for certain subsets of the models, e.g., the a-model, /3-model and 7/?y-models; one can also 
show this is the case for all models at very small and very large radius. Some models with 
a black hole have > (Tremaine et al. 1994) and the stability of these models is as 

yet unknown. 



3.2. Asymptotic expressions at small and large radius 



Asymptotic expressions are available for all (a, /3, 7)-models. In particular, the 
asymptotic expression for the surface light is also a double power-law, as for the volume 
density, but with a slope 1 — /3 at large radii and 1 — 7 at small radii (if 7 > 1); if 
< 7 < 1, the cusp only appears in the volume density. So the asymptotic behavior in the 
integrated light of the (a, /3, 7)-models are nearly the same as the set of double power-law 
models used by Kormendy (1994) to fit the observed galactic nuclei. The properties of the 
(a, /3, 7)-models are also qualitatively the same as the 7/?y-model at small radii and can be 
classified into three types as in Tremaine et al. (1994). But the models are much more 
general than the 7/?y-model at large radii, because the asymptotic volume density profile 
is rather than the fixed r~'^, hence they can fit a wider range of observations. The 
asymptotic expressions for cr^(r), crp[R) and I{R) are given as follows. A black hole of mass 
1^BH is included. 

At small radii, 



m 



(1 + 7)? 



if 7 < 1, 



, , 1 mBH 
Awr log - + — — — 
r (1 + 

rriBH 



if 7 = 1, 



+ 



(1 + 7)^ 
A3 if 7 < 1, 

An log— if 7 = 1, 

A^R^-^ if 7 > 1, 
rriBHAs 



if 7 > 1, 



As + 



Ai2_ rnBuAu 

1 "T 



log 



R 



i?logl 



if 7 < 1, 

if 7 = 1, 



(28) 
(29) 

(30) 
(31) 
(32) 
(33) 
(34) 

(35) 



R 
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A,R'-. + if 2 > 7 > 1, (36) 

-4,3ii4logi + ^ itT = | (37) 
A.if-M^ if7>| (38) 



where the constants Ai — A14 are given in the Appendix. 
At large radii, 

ni?) - (40) 

- ^^^^^ (41) 
where the constants D4 and Dg equal to A4 and Ag respectively if replacing 7 with /3. 

4. Summary 

In conclusion, I have presented a general family of models with analytical potential- 
density pairs which have asymptotic behaviors the same as observed galactic nuclei. More 
specificly, the density distributions have a double power-law radial profiles, which can be 
used to fit the central cusp, the power-law tail and the transition regions of observed galactic 
nuclei. In the spherical case the models generalize the 7/?y-models by Dehnen (1993) 
and Tremaine et al. (1994), and greatly enlarge the set of models with good analytical 
dynamical properties. Analytical expressions and asymptotic expressions for the intrinsic 
and projected velocity dispersion of the models are given for cases with and without a 
central black hole. 

The structure and dynamics of a general galactic nuclei are often more complex than 
in a spherical isotropic double power-law model. Both the non-spherical distribution in the 
light and a possible anisotropy in the velocity need to be modelled numerically, e.g., with 
the Schwarzschild (1979, 1982) approach and with the N-body approach. However, the 
analytical potentials are still very useful as they can be generalized to the non-spherical 
case and applied in these numerical methods as a fast solver for Poisson's equation. 



I thank R. Michael Rich and David N. Spergel for encouragements and comments on an 
earlier draft, Hans- Walter Rix, Simon White and Walter Dehnen for discussions and critical 
readings of the manuscript, John Hibbard for giving me the reference to the 7/?y-models. 
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5. Appendix 

In the following, I give detailed expressions for a few quantities used in the main text. 

1. The three parameters p;, qi and q in equation 7 and 11. 

pi = ai{--fi -1 + 2) ,qi = ai{/3i -1-3) ,ci = ai{/3i - -fi) . (42) 

2. The incomplete Beta- function i?(a, 6, x). 

B{a,h,x)= r dtr-^{l -tf-^ . (43) 



It is a simple integration, which reduces to elementary functions of x if a or 6 or 
1 — a — 6 is a natural number. One can prove this for the case that 6 is a natural number by 
expanding the polynomial (1 — t)^~^ and performing the integration for each term. Similarly 
for the other two cases after a change of variable from ttol — tor(l/t — 1). 

In cases that the function cannot be reduced to elementary ones or the expression is 
too lengthy, it can be numerically computed by a function call to, e.g., the efficient BETAI 
program in Numerical Recipes (Press et al. 1992); it converges in 0[^J Max[a^ b)) iterations. 

3. The normalization constant C . 

C = i . (44) 

47ra5(a(3 -7),a(/3 - 3),1) ^ ' 

4. The constants Kn,h ^n,i and w in equations 12 and 13. 



4(n + w)^ - 1 
IGvra 

1 2*'^+^a{n + w) nlV^iw) 



Kn,l = ^ 7._.2 ' (46) 



(47) 



Nn,i Tr{4:{n + w)^-l)T{2w + n) 



where 



w = {2l + l)a + ^. (48) 



5. The constants Ai — A14 in equations 28 to 38. 

A, = 47ra'C t x-'^-^^+^'^-^l - x)"^^+'^-'dx f ^^('-"'^-'(l - xi)"^'^-^)-^^! 
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_ 27rC 

(3-7)(7-l)' 

A3 = 2aCB{a{l - 7), - 1), 1), A4 = CB^^, ^ 1), 

i?(a(l-7),«(/3-l),l) ' 

^^-(3-7)i^(-^,il)' ^ ^ 



A9A4 , 5(Mii) . 1 

y-l 1 
2 ' 2' 



Aio , A3A5 , 8 A10 ^ A9A4 

1 1 1 \ ' 



- ^' - - 35(1,1,1)' "^^^ - ^ 



6. Si{x) in equations 14 to 21. 

S,{x) = -logx if^ = (50) 
I - x' 

= ; — otherwise. (51) 

t 

where x is given in equation 5. The function is singular at x = for z < 0. 

7. The coefficients a^, 6i, di — hi in equations 14 to 21. 

iiraC 



a(3-7) + z ' ' ^.^0 



q{a{/3 - 3) - b, = aJ2<l{(^ - ^J)(^^-J, d, = ^eja,_j, 

1=0 



= aCq{a{/3 + 1) - f^ = 2Cq{a{2--f),i), g, = ^hja,_j, = 2Cq{a/3,i), 

j=o 

■«!... 

q{i,i) = { — ^y~ri 77, if ^ ^ i ^ , otherwise = 0. (52) 

j!(« - Jy- 

8. The functions u\"\R) and Vi"\R) in equations 20 and 21. 

UI"\R) = r\ ^,,rl^"l. yd-, (53) 
"'0 + sm ' u 

= / (-— --^7— )'cos^^^J.. (54) 
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One can further reduce u\"\R) and V^"\R) to elementary functions for all values of i by 
the following recursion relation for the a = 1 sequence. 



anc 



where 



u\'MR) = u\'\R) + -^^^^j^, (55) 



II) + 

V!^{R) = v!'\R)-'i^^^^^-^, (56) 



I dR 
R dVl^\K) 
dR 



ui'\R) = -RZ{R) + {^-l)R, (57) 

vI'\r) = -{l-R^)Z{R) + {^ + l)R' -R-1, (58) 

V^'\K) = J, (59) 

V^iB) = ^ + 1, (60) 



Z{R) = RX{R) - 1 (62) 

X{R) = 1 if i? = 1, (63) 

cosh-Vl/i?) , , , 

arccos(l/i?) 

= (i^2 _ 1)1/2 Otherwise. (65) 



Note that near i? = 1, stable numerical evaluation of X[R) needs to make use of the the 
Taylor expansion given in Hernquist (1990). 



Table 1: Analytical models 


in the (a, /3, 7) -family 




Analytical quantities 


parameters (a,/3,7) 


Comments 




(a,/3,7)* 


non-spherical: eq. 8 and 9 


$(r, 6, tb) 


(n,3 + -,7) 


eq. 17 




and others^ 


non-spherical: eq. 10 and 11 


a^(r) and M(r) 


(n,3 + ^,7) 


eq. 19 and 16 with optional BH 


al{R) and I{R) 


(1,3 + A;,0) 
(1,3 + A;,l) 
(1,3 + A;,2) 


eq. 20 and 21 with optional BH 


HE) 


7/?y-models with 7 = 2 — - 


Dehnen 1993/Tremaine et al. 1994 




(1,5,0) 


Plummer (1911) 




(1,5,1) 


eq. 26 without BH 




(2 ^ ^) 


eq. 27 without BH 


Special cases: 






a-models 


K3 + ^,2-l) 




/3-models 


(1,/3,1) 


«('■) = -i-i+?^ 


7/?y-models 


(1,4,7) 


Hernquist model (1990) 


(1,4,1) 


cusped 


JafFe model (1983) 


(1,4,2) 


cusped 


Plummer model (1911) 


(1,5,0) 


finite core 


Perfect sphere 


(1,4,0) 


finite core (de Zeeuw 1985a, b) 


Mod. Hubble profile 


(1,3,0) 


finite core (Binney & Tremaine 1987) 


Mod. isothermal sphere 


(1,2,0) 


finite core (Sackett and Sparke 1990) 



* a, /3 and 7 are the transition width at the break radius r = 1, the slope of the outer and 
inner power-law respectively. Also G = M = 1. 



^ If n and k are any natural numbers, 1, 2, 3, etc, or if one of the four pairs (q — pi, qi), (pi, 
Q — <?/), {Ph <?/) s-iid (q — Pi, ci — qi) is a pair of natural numbers; p;, qi and q are defined in 
eq. 42 and / is the angular quantum number. 



